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. In this letter we discuss cosmological first order phase transitions with de Sitter bubbles nucleating on 

Q ' (inhomogeneous) matter backgrounds. The de Sitter bubble can be a toy model for an inflationary phase 

y-^- ' of universes like our own. Using the thin wall approximation and the Israel junction method we trace the 

, classical evolution of the formed bubbles within a compound model. We first address homogeneous ambient 

■ space (FRW model) and already find that bubbles nucleated in a dust dominated background cannot expand. 

For an inhomogeneous dust background (LTB model) we describe cases with at least initially expanding 
bubbles. Yet, an ensuing passage of the bubble wall through ambient curvature inhomogeneities remains 
unnoticed for observers inside the bubble. Notable effects also for interior observers are found in the case 
of a rapid background phase transition in a FRW model. 
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1 Introduction 

> , 

The influence of background inhomogeneity can for instance be important in the context of string landscape 
sampling via tunneling processes. Standardly, bubble evolution has been studied almost exclusively in de 
Sitter spacetime, see e.g. [1], given the long lifetimes of the metastable vacua in the landscape and the 
cosmic no hair conjecture. However, there are recent scenarios in which tunneling is catalysed and can be 
made very rapid, such as chain inflation [2] or DBI and resonance tunneling [3, 4]. Therefore we consider 
' bubble propagation on backgrounds that contain homogeneously and inhomogeneously distributed dust. 

. Inhomogeneous initial states could be relevant especially for processes of resonant tunneling, see [5, 6] 

and [7]. Another context is the question of whether and how non-standard backgrounds influence the onset 
of inflation, see e.g. [8]. More details on the methodology, more numerical examples and a perspective on 
the tunneling process on time dependent backgrounds can be found in [9]. 



2 Junction Method 

The bubble wall, which we assume to be thin w.r.t. the bubble size, separates interior and exterior space- 
times. In order to analyse the effects of a matter background on the bubble evolution we make use of 
the Israel junction conditions [10], similar to [11]. The interior spacetime is fixed to be de Sitter and can 
be a toy model for the inflation of our own universe. We describe the bubble wall by the induced met- 
ric hijdy'dyi = — dr 2 + R 2 d£l 2 and stress-energy = —ahij, with bubble surface tension a. The 
non-standard part is the exterior which is modelled by a Lemaitre-Tolman-Bondi (LTB) model [12]: 

,9 ,9 (rd r a(t,r) + a(t,r)) 2 , 9 . 9 , 

d.s 2 = -dt 2 + V 1 \ " dr 2 + a 2 (t, r r W . (1) 

1 + 2ki(r) 
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The LTB spacetime is spherically symmetric but radially inhomogeneous. The function E(r) is the local 
curvature of constant time slices. With dust plus vacuum source, two non-trivial field equations remain 

d t a\ 2 2E 2M A n 2d r M 

~ T2 = + ~q and 8n P = 2 21 a ! — T ' (2) 

a J a z r z a^r 3 3 a z r z [ra r a + a) 

There are three free functions in the model: M(r), E(r) and the bang time. At the same time there is 
a remaining gauge freedom in r. For d r M > one can rescale r such that M(r) = ^-Ar 3 , with A a 
constant. So, the LTB model leaves us with two free functions: the local curvature and the bang time. 

In order to glue the spacetimes together we use the method of Israel [10]. This implies two junction 
conditions for the discontinuity of the induced metric and extrinsic curvature on the bubble wall S: 

[hij] = h±\z - h^\z = and [Kij] - hy [K] = 87^ , (3) 

with ± labels denoting exterior/interior spacetime. Note that the second equation is closely related to a 
continuity equation of the form V^S"* + \T^n a e-] = 0, where n M is a normal vector on E. We consider 
the exterior stress-energy as given by a perfect fluid and the interior one by a cosmological constant. 
From the junction conditions we get the final equations of motion in the LTB frame, they read 



a - -(1 + 2E)fd t a + y/ 1 + 2E) 1 + 2V) Y (fdra + a)d t f 

dtr = 7=a 1 — VT^ — ^TT\ > dt<7 = P 1 (4) 

(rdra + a) {2E - 2V) , / , , ,2 m -,2 



1 + 2E- (fdfa + af (d t r, 
with r (t) the bubble trajectory in LTB coordinates and Y = (rdta) 2 — 2E + 2V. The potential V reads 

"A_ (A A + -A_ 



2V 



3 V 3a 3 o- 24na 



2no 



9 / 8ttA A, - A_ 
R 2 with Ana < ^/ — + -t- ( 5 ) 



being an additional and purely geometrical constraint for the surface tension that stems from our choice of 
fixing the sign of the normal vector on E. 

First, let us consider the homogeneous case. Already the study of this case can put restrictions on the 
applicability of the described junction method and will throw off some interesting physics. In addition 
to A± we put an amount of dust density in the outside space. Its initial distribution is assumed to be 
homogeneous and local curvature is assumed to vanish: this is the homogeneous limit of the LTB model 
which reduces to a FRW model. The fate of a bubble of new vacuum that nucleates comovingly (dtr = 0) 
on such a background will depend crucially on the force balance and thus on the relation of the latent heat 
to surface tension and dust density. After differentiating (4) we have 

„ 2 1 1 /A+ -A_ „ 2p\ 

d?f\ a n = - — 2-ko -) . (6) 

* a \ 247rcr 3cr J 

In contrast to nucleation on a pure vacuum background the presence of matter can now impede the ex- 
pansion of a vacuum bubble. There is a competition between the surface tension that tries to collapse the 
bubble and the pressure difference from the latent heat of the vacuum. In the pure vacuum case the bubble 
wall will move in the direction of smaller pressure, i.e. into the space of higher vacuum energy. As soon as 
a non-negligible contribution of pressureless dust is added to the budget the pressure support cannot prevail 
anymore and bubbles would airways collapse. Only for vacuum dominated setups (p < p vac = A + /8tt, 
and A + > A_ > 0) expanding bubbles can exist. In a dust dominated background the initial acceleration 
of the bubble becomes negative leading to bubble contraction. The situation is illustrated in Fig. 1. 

Now let us consider the inhomogeneous LTB case. There are two functions with which we can introduce 
inhomogeneity: the local curvature and the bang time. A choice of the latter is equivalent to a choice of 
the initial dust density profile po(r) = p(to,r). Thus we have E(r) and po(r) at our disposal. What we 
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Fig. 1 The early-time behaviour of comovingly nucleated de Sitter bubbles on FRW background with vacuum as well 
as dust source. The fate of the bubble as seen by an exterior comoving observer crucially depends on how the dust 
density compares to the surface tension of the bubble a and the latent heat of the vacuum e vac = (A+ — A_) /8n, 
see also eq. (6). In order to to produce an expanding bubble only little dust density in the background is allowed - on 
matter dominated backgrounds (p > p v ^ c — A+/87r) the nucleated bubbles must collapse. The shaded region in the 
plot is excluded by the junction constraint (5). 

are looking for in the end is whether the motion of a bubble through an inhomogeneity leads to potentially 
observable disturbances in the trajectory of the wall. The result of the last section shows that the presence 
of dust has significant influence on the trajectory and makes it difficult for the bubble wall to propagate 
to an exterior inhomogeneity at all. Note that also relaxing the condition of a comoving nucleation and 
considering bubbles with 9 t f(to) > does not save the bubbles from collapsing whenever there is enough 
dust in the background. However, we can enforce a transition by combining both free LTB functions. 
First we take an inhomogeneous initial dust distribution with a radially increasing profile - then small 
bubbles nucleate in a still vacuum dominated region and can initially expand. Second we take a curvature 
inhomogeneity and place it right in the way of the expanding bubble wall - note that the curvature profile is 
constrained due to shell crossing. For this twofold inhomogeneous setup we have compared the evolution 
of the bubble wall trajectory with its evolution for the homogeneous case and found a difference only in 
outside coordinates. In the (more relevant) inside frame the transition remains unnoticed, cf. Fig. 2. 

We also considered an FRW background with a perfect fluid that is just undergoing a phase transition, 
say a reheating or vice versa phase transition, while de Sitter bubbles evolve on it [9]. In contrast to the 
dust case the phase transition leaves a distinct perturbation of the wall trajectory, also seen in the inside. 

3 Conclusions 

On a de Sitter background vacuum bubbles generically expand into the space of higher vacuum energy. 
Here we have seen that a dust dominated homogeneous background hinders the expansion of de Sitter 
bubbles. Consequently, vacuum bubbles can hardly propagate towards given matter inhomogeneities in 
the background. In the LTB model we have used for the inhomogeneous background we have two free 
functions at disposal: the initial dust density profile and the spatial curvature profile. For a special choice 
of these functions the bubble wall reaches the inhomogeneities but these do not affect the trajectory of the 
wall but merely the surface tension, as seen by interior (physical) observers. Nevertheless, we found that 
also the interior observers can potentially see a significant perturbation of the bubble wall if the curvature 
inhomogeneity in the ambient LTB model is replaced by a smooth phase transition, say w = 1/3 — > w = 
— 1 or vice versa, taking place in a FRW background, cf. [9]. 
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Fig. 2 The numerical evolution of a growing de Sitter bubble on an LTB background which is inhomogeneous in its 
initial dust profile and has a curvature inhomogeneity. The profile of the latter is shown in the upper left figure: its 
form is chosen as a tanh, its amplitude is bounded by shell crossing. The upper right figure compares the trajectories 
of the vacuum bubble on an undisturbed flat background (red) with such on the LTB background (blue) as seen from 
the (exterior) LTB frame. The distinct perturbation from the curvature disturbance seen in this frame vanishes when 
the trajectory is computed in de Sitter coordinates (lower left figure). The effect is a coordinate artefact - nevertheless, 
considerable fine-tuning is needed in order to be even able to study such transitions within the junction method. As a 
consequence of the junction approach there is a notable effect on the surface tension of the bubble (lower right figure). 

If such disturbances of our bubble wall do exist it might not be impossible to observe them. In [ 1 3, 14] a 
different mechanism - primordial bubble collisions - has been discussed and it is argued that perturbations 
in the bubble trajectory can actually lead to a characteristic modification in the redshift of the reheating 
surface and thus could potentially be detectable in the CMB. 
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